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1 Introduction 



There has been a remarkable progress recently in understanding the two-dimensional 
quantum gravity. There are two main motivations to study the two-dimensional 
quantum gravity coupled to matter. Firstly, it is precisely a string theory when 
the two-dimensional space-time is regarded as the world sheet for the string. Sec- 
ondly, it provides a toy model for the quantum gravity in higher dimensions such as 
four dimensions. There are two approaches to study the two-dimensional quantum 
gravity. First one is the matrix model which gives a discretized version of the two- 
dimensional quantum gravity 0. The second one is the Liouville theory as 
a continuum theory 0, 

The matrix model is quite powerful in providing a nonperturbative treatment. 
Exact solutions of the matrix model have been obtained for two-dimensional 

quantum gravity coupled to minimal conformal matter with central charge c < 1. 
The result can also be understood by means of the continuum approach using the 
conformal field theory [Q, In both approaches, it has been very difficult to 
obtain results for two-dimensional quantum gravity coupled to conformal matter 
with central charge c > 1. In matrix model approach, it is easy to write down matrix 
models for cases with c > 1 as well However, these models are not solvable up 
to now [0. The numerical simulations suggest that it is not at all obvious if a 
matrix model candidate to describe a c > 1 model has the continuum description 
[p. Although several exact solutions of the matrix model have been obtained, it 
is worth studying approximation schemes which enable us to understand unsolved 
matrix models at least qualitatively, especially for the case of c > 1. In order to make 
use of such a scheme, however, we first need to make sure that the approximation 
method gives correct results for the exactly solved cases. 

Recently Brezin and Zinn- Justin have proposed a renormalization group ap- 
proach to the matrix modelQ. They drew an analogy between the matrix model 
and the critical phenomena in statistical mechanics. The matrix model possesses 
the double scaling limit 0. The limit specifies how the coupling constant should 
approach to a critical value as the size of the matrix N goes to infinity. This limit 
can be regarded as a continuum limit where bare parameters should approach to 
a critical value as the cut-off of the theory goes to infinity. The critical exponents 
determine how the coupling constant should be tuned to reach the continuum limit. 
Therefore the double scaling limit of the matrix model may be regarded as coming 
from the fixed point of the renormalization group flow. They observed that a change 
N —>■ N + SN can be compensated by a change of coupling constants g —>■ g + 6g in 
order to give the same continuum physics. They needed to enlarge the coupling con- 
stant space as in the Wilson's renormalization group approach Consequences of 
their approach have been examined by several groups |TT]. A similar approach has 
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been advocated previously for the expansion in a somewhat different context 



such as the 1/N corrections or the string field theory [12 . 

In the case of the one-matrix model with c < 1, Brezin and Zinn- Justin obtained 
reasonable results for the fixed point and the susceptibility exponents in the first 
nontrivial approximation. In order to demonstrate the validity of the renormaliza- 
tion group approach, however, one should show that the systematic improvement of 
their approximate evaluation converges to the correct result. To this end, it is im- 
portant to study a model in which a renormalization group equation can be derived 
exactly, even if it is simpler than the matrix model. The vector model has been 
proposed for a discretized one-dimensional quantum gravity, in the same way as the 
matrix model for a discretized two-dimensional quantum gravityp!3|-||16||. 

Recently we have analyzed the vector model by means of the renormalization 



group approach and have clarified its validity and meaning This paper is a 

more complete account of our analysis of the vector model and some results on the 
matrix model. For the vector model, we have obtained an exact difference equation 
which relates the free energy — \ogZ]si-2{g) to the free energy —\ogZ]si{g — 2Sg) 
with slightly different values of coupling constants. We have found that these cou- 
pling constant shifts 6gk are of order and occur in infinitely many coupling 
constants. We also obtained infinitely many identities which express the freedom 
to reparametrize the field space. By using these identities, we can rewrite the flow 
in the infinite dimensional coupling constant space as an effective flow in the space 
of finite number of coupling constants. The resulting effective beta function deter- 
mines the fixed points and the susceptibility exponents. The inhomogeneous term 
in the effective renormalization group equation serves to fix non-universal (analytic) 
terms of the free energy. To illustrate the procedure by an explicit example, we an- 
alyze in detail the cases of one and two coupling constants which give the first two 
multicritical points m = 2 and 3. We obtain the fixed points and the susceptibility 
exponents which are in complete agreement with the exact results. 

In Sect. 2, the renormalization group approach for large N models is described 
and our results on matrix models are given. In Sect. 3, the exact difference equation 
is obtained for vector models. In Sect. 4, reparametrization identities are obtained 
and the effective renormalization group equations are derived for the single coupling 
constant case. In Sect. 5, the case of the multi- coupling constant is worked out. The 
two coupling constant case is studied in detail, and the renormalization group flow 
is illustrated in the two-dimensional space of coupling constants. Sect. 6 is devoted 
to a discussion. 
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2 Renormalization Group Approach to Matrix Mod- 
els 



We first recall the renormalization group approach for the matrix model. The par- 
tition function Z^ig) of the matrix model with a single coupling constant g for the 
cubic interaction is defined by an integral over an x hermitian matrix $ 



ZN{g) = J rf^'$exp [-ATtr + 1$^) 



(1) 



Here we consider the cubic interaction, while Brezin and Zinn- Justin studied the 
case of the quartic one. It seems that the Z2-symmetry of the quartic potential is 
not particularly useful to perform the perturbation in higher orders. We take the 
cubic one because we need less calculation in that case. The matrix model gives 
the random triangulation of two-dimensional surfaces [0]. The expansion of the 
free energy F{N,g) 

-| oo 

m,g) = -^iogZ^(^?) = j:N-''Mg) (2) 

h=0 

distinguishes the contributions fh from the surface with h handles. In the double 
scaling limit 

2 

— ^ OO, g ^ g*, with A^Ti {g — g^) fixed, (3) 

the singular part of the free energy satisfies the scaling law ^ with the susceptibility 
exponent 70 -|- /i7i linear in the number of handles h 

fh{g)sing = {g-g*y-''''-'''''a, + ---, 

F{N,g),,^g = (^-^,)2-7o/(ivi(^_^,)). (4) 

The functional form of / can be determined by a nonlinear differential equation 
(string equation) 0. The continuum limit is achieved at the critical point g ^ g* 

where the average number of triangles diverges. The double scaling limit suggests 

2 

that we may draw an analogy between A^^i and the momentum cut-off 

g - g* ^ ^ ntriangics 00 
2 

A^Ti ^ 00 <^ a (5) 
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where the lattice spacing of the random surface is denoted as a and -v^ indicates 
the conjugate relation. There is also an analogy between the scaling variables to be 
fixed in the double scaling limit and the renormalized variables in the continuum 
limit 

2 

N''^ {g -g*) ^ a- ntriangles- (6) 

It has been proposed that the free energy F{N,g) of the matrix model satisfies 
the following renormalization group equation M 



F(N,3) = r{g), (7) 



where /3{g) is called the beta function. The anomalous dimension and the inhomo- 
geneous term are denoted as 7(5') and r{g) respectively. A fixed point g^ is given by 
a zero of the beta function. The exponents of the double scaling limit can be given 
by the derivative of the beta function at the fixed point 

To obtain the beta function, we integrate over a part of degrees of freedom of 
the (iV + 1) X (A^+ 1) matrix. 0, 0. We parametrize the (A^+ 1) x (A^ + 1) matrix 

^N+i as 

-(!;:). (10) 

where f is a complex A^-vector and a is a real variable. Then we integrate over v in 
the first order of the perturbative expansion to obtain 

Z^^,{g) = {XM{g)f"Z^{g + 6g), (11) 

where 

Sg = ^^^^^ + 0{N-% (12) 

A = l + ^^^ + 0(iV-^). (13) 
We can read off /3, 7 and r in (^ as 

m = ^(11^?'-^?), (14) 

l{g) = 2, (15) 
rig) = ^(3/-l). (16) 
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Here we have used the assumption that there is no contribution to these functions 
from the integration over a, which needs to be reexamined for our cubic potentials. 

We find a zero of the beta function at = — 1/vTl = —0.302 . . . and identify 
this with the critical coupling constant. Using 
as 

7o = 0, 7i = 2. 

Compared to the exact result 
= -(432)-!/^ 
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the critical exponents are given 

(17) 



-0.219..., 70 = -1/2, 71 = 5/2, 



this approximation seems to be reasonable. This situation is similar to the case of 
the quartic potential studied by Brezin and Zinn- Justin P|. In fact, the value of the 
critical exponent at this first nontrivial order is precisely identical to their result. 

In the first order approximation above, induced interactions tr $^ and tr <l>^ of 
order g'^ and g^ has been ignored consistently. To perform the next order calculation, 
we need to enlarge the coupling constant space and start with the partition function 



1 



-Nil { , 

V2 3 



h 
< 

4 



(19) 



Actually it turns out that, if we assume that f g'^,h g^, it is sufficient to include 
these two coupling constants for the consistent second order calculation. We obtain 
the beta functions 



f3g{9,h, /) 
Phi9,h, f) 
Pf{9,hJ) 

We find a common zero at 



-h-g\ 



ng^ + 18gh-6f), 



9* 



-3(3/ -2^?^). 



-9> 
2 



9I 



44 

-0.00617.. 
-0.00115 



-0.280, 



(20) 
(21) 
(22) 



(23) 
(24) 
(25) 



The scaling matrix at the fixed point 



/ d(3g/dg d(3Jdh d(3gldf \ 

df3uldg df3h/dh dPh/df 
\dPf/dg d(3j/dh dPf/df J 



(26) 



g=g*,h.=h,,f=f^ 
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exhibits three eigenvalues 1, —1.10 . . . , —1.55 . . .. Substituting the positive one into 
d^), we obtain the critical exponents 

70 = 0, 71 = 2. (27) 

Though approximate value for the critical coupling g^, improves slightly, those for 
the critical exponents show no improvement. Before performing more elaborate 
calculation for the matrix model, we find it more illuminating to study the case of 
the vector model where we can clarify the situation more fully. 



Renormalization Group Approach to Vector Mod- 
els 



The partition function of the 0{N) symmetric vector model is given by 

ZN{g) = 



d 



exp 



,2\k 



k=l 



(28) 



where is an dimensional real vector |T3[ - [|I6|. Here we introduce infinitely 
many coupling constants Qk, since we need all possible induced interactions after a 
renormalization group transformation even if we start with a few coupling constants 
only. The 1/A^ expansion of the logarithm of the partition function gives contribu- 
tions from h loops as terms with iV^~'^. The vector model has the double scaling 
limit N ^ oo with N^^'^'^ {g — g^) fixed, where the singular part of the free energy 
satisfies the scaling law 



- log 



Zn{9) 



ZM{gi = l.gk = ^{k>2)) 



sing 



h=0 



g.)^-^'^-^^W + 



(29) 



One should note that the power of 1/A^ and the definition of the susceptibility expo- 
nents are slightly different from the matrix model. Therefore the relation between 
these susceptibility exponents and the derivative of the beta function becomes 



1 



71 



70 = 2- 



7(9*) 



(30) 



In the spirit of the approximation method of ref. 
(A^ + l)-th component a of the vector (pN+i 

0Af+l = (0Af,a), 
Zn+M = / d''<PN da exp -{N + 1) ^ ^' 



we can integrate over the 



k=l 



2k 



N 



+ a 



2\k 



(31) 
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Neglecting higher order terms in we obtain 



ZN+i{gi = l,gk = 0ik>2)) 



J (i^0Arexp 



gk + ^gk 



k=i 



2k 



where the shifts Sg^ of the couphng constants are found to be 



k=l 



k 



\k=l 



9k^k-i 
91 



k=l ^ 



+ 



(32) 



(33) 



Now we shall show that the shift 5gk can be evaluated exactly in the vector model 
without using approximation methods such as described above. To demonstrate that 
the result is exact, we start with the partition function Ziq_2{g) . After integrating 
over angular coordinates in R^~^, we perform a partial integration in the radial 
coordinate x = d)^ 



(g) 



TT 2 



■N 



IT 2 



iWo 



dx 



exp 



-(iv-2)E 



r(f- 



1) Jo u=l 



gkx 



„k-l 



k=l 



exp 



2k 



n 



.(34) 



Identifying the right hand side with {{N — 2)gi/27T)ZN{g — 26g), we obtain an exact 
difference equation for the logarithm of the partition function 



[(-logZ^ig)) - i-\ogZN-2ig))] - log 



{N-2)g, 



2tt 



-[(-\ogZ^{g-26g)) - {-logZ^)]- 



(35) 



We find that the shifts Sgk of the coupling constants are exactly identical to the result 
(|33|) of the approximate evaluation. We would like to stress that no approximation 
is employed to obtain eq. (p5|). Therefore we can infer that the vector model offers 
an example to justify the approximate evaluation method of ref.0. 

In the N oo limit, we can obtain a differential equation from the exact 
difference equation (^) 



J;^(-logZ^(^))-^log^ 



oo g 

fc=i 99k 



log ZN{g)). 



(36) 



One can bring the quadratic term in the potential to the standard form 0^/2 since 
gi can be absorbed by a rescaling gicf)"^ 0^. We have 



ZN{gug2, g3, ■■■) = gi ^^^^7v(i, g2/gl, gs/g't, ■■■)■ 



(37) 
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Therefore it is convenient to use the rescaled couphng constants gk together with gi 
as independent couphng constants 



9k = gk/gt 

We shah define the free energy for the vector model 

F{N,~g) 



1 1 1 , Ngi 
logZATigf) log , 



{31 



(39) 



If we use gi,g2,g3,---a.s independent coupling constants, we find from the rescaling 
identity ( P7D that the partition function Zj\f depends on gi only through the factor 
gi ^^"^ . Therefore, the free energy F{N,g) is independent of gi and is a function of 
h {k > 2) only. 

We denote the partial derivatives with respect to (71,^2,^3, ■ ■ ■ by |g, and those 
with 5-1, 5-2, 5-3, ■ ■ ■ by Ig 



d_ 
dgi 



_d_ 

dgi 



g k=2 i 



d 




d 




1 d 


dgk 


9 


dgk 


9 


9i dgk 



Thus we obtain a renormalization group equation for the free energy F 



N 



_d_ 

dN 



fc=2 \9i 9i J 09k 



+ 1 



F{N,~g) = N5g^ 



2gi 



Eq.(41) shows that the anomalous dimension is given by 

1C9) = 1, 

which implies a relation between two susceptibility exponents 

7o + 71 = 2. 

We read off the beta functions in the rescaled coupling constants as 

m) = A' f % - ^**) . 

V 9i 9i J 

These beta functions can be evaluated explicitly using eq . (p3|) as 

f32{9) = -^2-3^2 + 2^3, 

Psig) = -2^3 + ^2 - ^hh + 3^4, 



f34{g) = -3^4 -92+ 4^2^3 - 2^3 - 8^2^4 + 



5; 



(40) 



(41) 



(42) 



(43) 



(44) 



(45) 
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Let us investigate the simultaneous zero of the beta functions /?a: = (A; > 2) in 
the spirit of ref.|0. The first condition (^2 = gives ^3 in terms of g2- The second 
condition (3^ = determines ^4 in terms of ^2 and ^3, and hence in terms of g2- As 
can be seen in eg. (^S)) , Pk{g) turns out to be a sum of kg^^i and a polynomial in 
gj,j < k. Therefore we find that there always exists a solution of = (fc > 2) for 
each given value of the coupling constant g2- In the following, we shall show that 
this strange result of the apparent existence of the one-parameter family of fixed 
points is due to a misinterpretation of the renormalization group flow. 



4 Reparametrization Identities 



The key observation is the ambiguity to identify the renormalization group flow in 
the coupling constant space. Though the above equation ( ^T]) seems to describe 
a renormalization group flow in the infinite dimensional coupling constant space, 
the direction of the flow is in fact ambiguous because all the differential operators 
[d/dgk) are not linearly independent as we will see shortly. To see this, note that 
the partition function ( PBD is invariant under reparametrizations of the integration 
variable (p. Since the model is 0{N) invariant, we can obtain new informations 
only from reparametrizations of the radial coordinate x = 0^. Since the radial 
coordinate should take values on the half real line, we consider the most general 
reparametrization which keeps the integration range [0, 00) 



X 



j=0 



(46) 



where e/s are infinitesimal parameters. Substituting (^) in ( plD and differentiating 
with respect to e/s, we obtain a family of identities 



LjZN{g) = 



=j+i 



U > 0), 



(47) 
(48) 



The differential operators Lj constitute half of the Virasoro algebra. We can show 



that this algebra is identical with the one found in |]T3| and ||T4 



It is more useful to use the rescaled coupling constants cjk- The reparametrization 
identity corresponding to eo is nothing but the infinitesimal form of the rescaling 
identity (^) and reads 

^ F{N,~g) = 0. (49) 
9 



dgi 
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In terms of the rescaled coupling constants, the reparametrization identity corre- 
sponding to £i is given as 



The reparametrization identities corresponding to is given by 



FiN,~9)=0. (50) 

9 



2j\ . d ~ . 




F(N,3) = 0. (51) 



We see that derivatives of the free energy in terms of infinitely many coupling con- 
stants Qk are related by infinitely many reparametrization identities. Thus one can 
expect that only a finite number of derivatives are linearly independent. The exact 
difference equation (|S3| ) combined with the reparametrization identities (^) and 
(|51|) constitute the complete set of equations to characterize the renormalization 
group flow in our approach. 

To illustrate the use of the reparametrization identities, we shall first take the 
case of a single coupling constant. Let us consider a point in the coupling constant 
space 

(5-1, h, 93, 9i, ■■■) = (91,92, 0, 0, . . .). (52) 

At this point, the j-th identity relates dF/dgj^2 to dF/dgj+i except for the case 
of j = 1 where an extra constant term is present. Therefore we can express 
dF/dc/k [k > 3) in terms of dF/dg2 by solving these reparametrization identi- 
ties recursively in the one-dimensional subspace (0) of coupling constants. It is 
most convenient to organize the solution in powers of We find explicitly at 

leading order 

dF BF 

Bk^ + Rk + 0{N-'), (53) 



where 



dgk dg2 



B^ = ^4^ (54) 





A = 4^2 + 1. (56) 
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These solutions at leading order are sufficient to eliminate dF/dgk {k > 3) in 



favor of dF/dg2 in (41). Thus we obtain a renormalization group equation with the 



effective beta function (3 {g2) and the inhomogeneous term r{g2 



where 



dN 



_d_ 

dg2 



/3^^(^2)— + 1 



F{N,~g2)=r{~g2), 



(57) 



/3^^(^2) 



rig2 



Bkf3k{g2, gk = 0,k>3) 

k=2 
1 

4 



A 



1-A£ /I 
1=^ log — 



2 1\ 32 1\2 

JSl I CO ~ ^ ^ 

-77^91 - 9 + Rkf3k{92, 9k = 0,k>3) 

^ ^ k=2 



(5^ 



1 



2VA 



1 1,1 1 i\ 2 / r, 

4 + 2^°^2 + 3(^^ + 4) + T5(^^ + 4' + 



(59) 



The graphs of the effective beta function P'^^{g2) and the inhomogeneous term 
r{g2) are shown in figure l.We find the beta function has two zeros g2 = —1/4 and 
92 = 0. 

The fixed point ^2 = ~l/4 is an infrared-unstable fixed point. By choosing the 
bare coupling constant around this point, we obtain a non-trivial theory with the 
renormalized coupling constant — cxd < ^2 < 0. Here we use the term "infrared" for 
the limit — > in analogy with the ordinary field theory. 

Furthermore, we can calculate the susceptibility exponents from the derivative 
of the effective beta function by using (|30D and (^) 



1 



7i 



f3'{g2.) 



3 
2' 



7o 



1(92* 

P'{92> 



1 

2' 



(60) 



The fixed point and the susceptibility exponent are in complete agreement with the 
exact results for the m = 2 critical point of the vector model corresponding to pure 



gravity 0-0. 
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Another fixed point ^2 = is an infrared-stable fixed point or the "gaussian" 
fixed point, since djS^^ /dg2 = — 1 < 0. Without particular tuning of the coupling 
constant we obtain a free theory corresponding to this fixed point. 

We observe that both the effective beta function /5'^^(^2) and the inhomogeneous 
term r(^2) are analytic in ^2 around the fixed point ^2 = —1/4. However, they 
become complex for ^2 > 0, as shown in fig.l. 

We can extract the complete information from the renormalization group flow 
in our approach, namely the exact difference equation and the reparametrization 
identities by a systematic expansion of the free energy in powers of 1/A^ 

00 

m^2) = E^"Vh(^2). (61) 

The free energy at leading order /o satisfies an ordinary differential equation exactly, 

h{92)-P''\92)^{92)=r{g2), (62) 
092 

where the effective beta function and the inhomogeneous term are those given in 
eqs.([58|) and (^Qf). We see immediately that the general solution is given by a 
sum of an arbitrary multiple of the solution of the homogeneous equation and a 
particular solution of the inhomogeneous equation. Since both the effective beta 
function ^'^^{92) and the inhomogeneous term r(^2) are analytic in ^2 around the 
fixed point ^2 = ~l/4, the singular behaviour of /o comes from the solution of the 
homogeneous equation. More explicitly, we can write /o as the sum of two terms 
fo{g2)sing and f 0(92) analytic- The former is singular around ^2 = -1/4 and satisfies 
the homogeneous equation, while the latter is analytic around ^2 = —1/4 and is a 
particular solution of the inhomogeneous equation (See figure 2 ). 



^ , , 1 , 1 + V A V A 

fo{92)s^ng = l^^&Y^'W^y ^^^^ 
f 0(92) analytic = ^ log + 2(1 ^ A) ~ ? ^^^^ 

It is important to notice that the singular term corresponding to the continuum 
physics (the so-called universal term) is specified by the beta function alone. Its 
normalization, however, cannot be obtained from the renormalization group equa- 
tion. The analytic contributions are determined by the inhomogeneous term and 
the effective beta function. 

The exact solution for /i = is obtained by choosing the normalization constant 
for the singular term f 0(92) sing to be unity. 



foih) = f 0(92) sing + fo{9: 



2) analytic} 
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1, 1 + VA 11 
= - log — + ^ - -. 65 

We find that the singular part of the free energy acquires an imaginary part for 
^2 > 0. this is due to the fact that the separation of the singular and analytic parts 
is defined by the behavior at g2 = -1/4. Both /o(^2)smg and foih) analytic become 
complex after extrapolation beyond ^2 > 0. However, the full free energy /o(^2) is 
of course real for ^2 > 0, in accordance with the stability of the potential in that 
parameter region. On the contrary, the full free energy /o(^2) acquires an imaginary 
part in ^2 < —1/4, which comes entirely from the singular part /o(^2)smg- This 
phenomenon corresponds to the instability of the potential in this parameter region. 

To obtain the full information for the free energy up to the order A^"'', we should 
write down the solution of the recursive reparametrization identities (0) and ( pID 
up to the order A^"'' explicitly. By inserting the solution into the exact difference 



equation ( p5| ) and by expanding it up to the power N , we obtain an ordinary 



differential equation for 



(1 - h)M~92) - P''T92)T^i~92) = rnih). (66) 

We see that the effective beta function is common to all h, whereas the inhomoge- 
neous terms Vh depend on h. Therefore we find the singular part [] of the free energy 
is determined by the beta function up to a normalization 

fhih) = fh{g2)sing + fhih) analytic^ 
fh{~92)s^ng = {~92 " 92*?"'''-'"'' + " " " • (67) 

If we sum over the contributions from various h, we find that the renormalization 
group equation determines the combinations of variables appropriate to define the 
double scaling limit. On the other hand, the functional form of the scaled variable 
is undetermined corresponding to the undetermined normalization factor at for the 
singular terms of each h 

oo oo 

F(iV,^) = Y: N-^U~9,)s,^, = Y: N-\~92 - ~92*f-''-''W 

h=0 h=0 

= {92 - ~92*?-''f{N'/^'i92 - 92.))s.ny (68) 

It has been known in the exact solution of the vector model that a nonlinear 
differential equation determines the functional form / [N^^'^'^(g2 — 92*)) of the 

V / sing 

^ We have included the possible logarithmic terms and integer power terms into the analytic 
contribution, since they are non-universal. The singular contributions are universal 
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scaled variables. The equation is called the string equation, or the L_i Virasoro 



constraint 0, [jT8[. The above consideration means that the complete set of our 
renormalization group flow equations, namely the exact difference equation and 
the reparametrization identities, does not give the L_i constraint for the universal 
singular terms. This is in accord with our objective of the renormalization group 
approach: to obtain at least fixed points and critical exponents even if the exact 
solution is not available. On the other hand, our equations give informations on 
the non-universal analytic terms. We would like to point out, however, that we can 
determine the functional form of / (N^/"'^(g2 — ^2*) ) if "we know the boundary 

V / sing 

condition of the renormalization group equation. Namely, we have only to know the 
functional form of F{N = Nq, ^2), the free energy with a fixed A^, or F{N, ^2 = ^2); 
the free energy with a fixed coupling constant. 



5 Renormalization Group Flow in the Multi- 
Coupling Constants Space 

It is known that there is a series of critical points labeled by an integer m = 2, 3, 4, . . . 



in vector models|jT3|-|jT5|. The critical exponents of the m-th critical point are given 
by 

7o = 1 — 1/m, 7i = 1 + 1/m. (69) 

It is known that the m-th multicritical points constitute a submanifold with codi- 
mension m — 1 in the infinite-dimensional coupling constants space. In other words, 
we have to tune m — 1 parameters to achieve the m-th multicritical point. 

In the previous analysis, we have found only the fixed point corresponding to 
the m = 2 critical point. It is natural that we have not found higher multicritical 
points, since we have limited our consideration in the one-dimensional coupling 
constant space. We can extend the analysis to the more general situation of finitely 
many coupling constants. Let us take 



9i 



1, 9k 7^0 {2<k<m), gk = {k>m + l). (70) 



In this subspace, each reparametrization identity involves only finite number of 
derivatives. Moreover, the j-th identity relates derivatives in terms of coupling con- 
stants ^j+i, ^j+2, ■ ■ ■ , gj+m- Therefore we can solve this identity to rewrite dF/dgj+m 
in terms of dF/dgk (j + 1 < ^ < j + m — 1). By successively using these identi- 
ties, we find that there are precisely the necessary number of recursion relations to 
express dF/dgk (A; > m + 1) in terms of dF/dcjk (2 < k < m). Consequently, we 
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can reduce the renormalization group equation effectively in tlie space of tlie finite 
number of coupling constants ^fc (2 < /c < m) 



F{N, 92,..., Qm) = r{g2, ...,gm). (71) 



Tfie m-tfi multicritical point sfiould be obtained as a simultaneous zero of all the 
beta functions [3^ = . . . = = 0. The susceptibility exponent is given by an 
eigenvalue of the matrix of derivatives of beta functions 

= dpf/d~gj (72) 

at the fixed point, which is an (m — 1) x (m — 1) real matrix. 

To illustrate the multi-coupling case explicitly, we take the case of two coupling 
constants. We can solve the reparametrization identities recursively to leading order 
in 1/N, and combine Pk{g2,93,gk — ^ {k > A)) to obtain the effective beta functions 

oeS(~ ~ N ^ -2a/37 {(3^ - 7^) log a + cychc a + + 7 /„on 
^^''^'^ (a/5 + /37 + 7«)M«-/5)(/3-7)(7-«) «/5 + /37 + 7«' 

r,eS(~ ~ ^ ^ 3a/j7 (/j^ -7^) log Q + cyclic 2 

^^''^'^ («/3 + /37 + 7«)M«-/3)(/3-7)(7-a) «/3 + /37 + 7« 

where a, (3,j are the three roots of the cubic equation 

Qsx^ + g2X^ + X - 1 = 0. (75) 

We find three fixed points as the simultaneous zeros of and (See figure 

3). 

The first fixed point is at (^2,^3) = (~l/3?l/27). This result agrees with the 
exact value of the coupling constants at the m = 3 critical point. We have evalu- 
ated the derivative matrix Qij = djdf^ / dgj by expanding the effective beta function 
around the fixed point 

a«(*^-is,^i)^(l I). (76) 

One of the eigenvalues of the derivative matrix is 3/4 which gives the exact suscep- 
tibility exponent 71 = 4/3, while the other eigenvalue is 1/2 which gives an analytic 
term. The corresponding eigenvectors also agree with the exact solution. As we 
stressed before, we can regard N^/'^'^ as an ultraviolet cut-off in taking the double 
scaling limit. Since both eigenvalues are positive, this is an ultraviolet fixed point. 
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or twice unstable fixed point. This means that we have to tune two bare couphng 
constants for taking the double scaling limit. 

The second fixed point is at ^2 = —1/4 and = 0. We obtain the derivative 
matrix Qij = dpf/dgj at (^2,^3) = (-i,0) 

n,, (92 = = 0) = (^l A), (77) 

which has eigenvalues 2/3 and —1/2. The eigenvalue 2/3 correctly gives the m = 
2 susceptibility exponent 71 = 3/2. The fixed point is infrared- repulsive in the 
direction {6g2, Sg^) oc (1, 0) and infrared-attractive in the direction {Sg2, Sg^) oc 
(—2, 1), which are the eigenvectors for the eigenvalues 2/3 and —1/2, respectively. 
The infrared-attractive direction is tangential to the m = 2 critical line 

~93 = -^[2 + 9^2 + 2(1 + 3^2)'/'], (78) 

~g3 = -^[2 + 9^2-2(1 + 3^2)'/'] (^2<0). (79) 

Moreover, a branch of the m = 2 critical line ([78| ) exhibits to be a trajectory of 
the renormalization group flow approaching to (0, —1/4) in the infrared limit. This 
means that we have only to tune the bare coupling constants on this line to achieve 
the m = 2 critical point. The renormalized coupling constants can take values 
between (—00, 0) and (0, 0) on the ^2 axis. 

The third fixed point is the trivial fixed point at the origin (^2,^3) = (0,0). 

We admit that we do not understand the reason why tuning on the other branch 
([79|) leads us to m = 2 criticality. Though the line (^) is a trajectory of the 
renormalization group flow connecting two fixed points, the flow reaches the gaussian 
fixed point in the infrared limit. This might be due to the fact that the scaling matrix 
at the origin has a divergence of log g^ order. 

We note that the m = 2 fixed point and the trivial fixed point are in corre- 
spondence with the ones found in the analysis of the single coupling constant case. 
Actually we find that the ^2-axis is an invariant surface and the effective beta func- 
tion i32^{g2,0) on it is identical to the effective beta function P2^{g2) in eg. (p^) 
obtained in the single coupling constant case. This is not surprising. Two opera- 
tions, setting ^3 to be zero and solving the reparametrization identities recursively, 
are commutative. 

Now we consider the general case of the flow in the m — 1 dimensional coupling 
constant space = {(^2, ^3, • • • , ^m, 0, 0, . . .)}. This space includes the subspace 
rk{2 < k < m). We find that the effective beta functions for the additional coupling 
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constants vanish everywhere on F^, 



/3Si = 0,--.,/3^^ = 0, if ^,+i = ... = ^^ = 0. (80) 

The effective rcnormaUzation group equation with g^+i = ■ ■ ■ = g-m = ^ turns out to 
be identical with that obtained in the k — Idimensional case. Therefore any fixed 
point in is a fixed point in F^ \im> k. 

The behavior of the free energy near a fixed point is determined by the derivative 
matrix Q at that point. We denote the j-th eigenvalue of as \j and its eigenvector 
as Vji 

E Vjk^ki = E V^>^^(~9*) = (81) 

k=2 k=2 

The renormalization group equation tells us that the singular part of the free energy 
is given by 



/i=0 



NF,,^g{N,~9) = Y.'^h Nj2iJlVjk{gk-~9k*)y 



l-h 



(82) 

j=2 k=2 

To take a double scaling limit, fine tuning of the coupling constants is necessary. 
If we start near the fixed point and let N decrease (infrared direction), only those 
coupling constants corresponding to positive eigenvalues of Q can survive. Eventu- 
ally, the combination of coupling constants corresponding to the maximal eigenvalue 
dominates in this limit and the susceptibility exponent 71 is given by the inverse of 
this maximal eigenvalue. 

Correlation functions among the operators 7V(0^)^, . . . , N^cf^Y can be dealt with 
in Fj„ if /c < m. These correlation functions are given by the derivatives of the free 
energy F with respect to the coupling constants §2, . . . ,gk- 

'^Wr ■ ■ ■ ^{<Pr") =J---- J-NF{N,g). (83) 



2ki 2kn I c ^Qkx dg 

In taking the double scaling limit, we can tune coupling constants by fixing all the 
scaling variables 

N^[Vjk{~gk-h*)]^: j = 2,---,m (84) 

k 

corresponding to the j-th eigenvalue Xj of the derivative matrix Q. Then the oper- 
ators corresponding to these combinations of coupling constants obey a scaling law 
near a fixed point determined by the derivative matrix fl. The scaling behavior of 
a single point function is 

{tyk-/^i'i>r)-N'^, (85) 
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where Aj is its eigenvalue. The operator belonging to a positive (negative) eigenvalue 
is a relevant (irrelevant) operator and is amplified (deamplified) as — oo. 

The exact result shows that the m-th multicritical point exists in and is 
isolated from for k < m. Fine tuning of the m — 1 coupling constants is needed 
for the double scaling limit at this critical point. The knowledge of the exact solution 
in the vector model |T^-|T^ helps us to conjecture the nature of the renormalization 
group fiow in the m—1 dimensional coupling constant space. When we consider space 
Fm of m — 1 coupling constants, we should obtain the critical points by demanding 

/5f = --- = /5™ =0. (86) 

The k-th multicritical point for m > k > 2 can be obtained as a fixed point in this 
m — 1 coupling constant space F^- For the k-th multicritical point, the fine tuning 
of /c — 1 parameters is necessary to achieve the double scaling limit. Corresponding 
to these tuning parameters, there are k — 1 relevant operators which should be given 
by linear combinations of (0^)^, . . . , {(f^Y- All other operators (0^)'^"'"^, . . . , (0^)™ 
should be irrelevant near the critical point. In other words, the /c-th multicritical 
point in Tm{fn > k) has k — 1 repulsive directions and m — k attractive directions. 

6 Discussions 

We have applied the renormalization group approach for large proposed by Brezin 
and Zinn- Justin to vector models. We have found that we have to consider an 
effective renormalization group fiow in the finite-dimensional coupling constants 
space in order to reproduce the exact critical coupling constants and the exact 
critical exponents. In the procedure, the reparametrization identities have played 
a crucial role. We have obtained explicit expressions of beta functions in the case 
of a single coupling constant and two coupling constants. We have found that the 
expressions obtained by approximate evaluation, namely the 1/A^-expansion and 
the perturbative expansion in the coupling constant, gives the result in complete 
agreement with those obtained from the exact difference equations. 

The phase structure for two-term potential has been clearly understood in the 
language of the renormalization group. As we described earlier, a naive application 
of the renormalization group approach to matrix models has given us results which 
show no improvement as we go to higher orders. In view of our analysis of the exact 
result on the flow of the coupling constants for the vector model, one should not 
be surprised by these results on matrix models. We should look for an effective 
renormalization group flow by taking account of the reparametrization freedom in 
the fleld space. Although the above results are derived exactly only for the vector 
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model, we have a similar ambiguity to identify the renormalization group flow in the 
coupling constant space for matrix models. For instance, there are identities which 
form a representation of Virasoro algebra in matrix models |T^. There has been 



some work to interpret the Virasoro constraints as reparametrization identities in 



the field space in the matrix model [|T9| . We are trying to devise a way to implement 
our idea of the effective renormalization group flow through the application of the 
reparametrization identities to matrix models in order to make the renormalization 
group approach more useful. Work along this hne is in progress. 

We thank B. Durhuus and T. Hara for an illuminating discussion and P. Crehan 
for a careful reading of a part of the manuscript. This work is supported in part by 
Grant-in-Aid for Scientific Research (S.H.) and Grant-in-Aid for Scientific Research 
for Priority Areas (No. 04245211) (N.S.) from the Ministry of Education, Science 
and Culture. 
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Figure captions 

figure 1 The effective beta function P'^^{g2) and tlie inliomogeneous term r{g2). 
Tlie real parts are sliown in solid curves and the imaginary parts in dashed curves. 

figure 2 The free energy in the leading order in N. The singular part and analytic 
part around ^2 = —1/4 are drawn separately. The real part is drawn in solid curves 
and the imaginary part in dashed curves. 

figure 3 Renormalization group flow in the infrared (decreasing A^) direction. One 
shows the global behavior of the flow, while the other shows the detail of the flow 
near the origin. The m — 2 critical line is shown in a soiled curve. 
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